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Update on EMD and Hilbert-Spectra Analysis of Time Series 

Goddard Space Flight Center, Greenbelt, Maryland 


U.S. Patent 6,381,559 presents fur- 
ther information about the method de- 
scribed in “Analyzing Time Series Using 
EMD and Hilbert Spectra” (GSC- 
13817), NASA Tech Briefs, Vol. 24, No. 10 
(October 2000), page 63. To recapitu- 
late: The method is especially well 
suited for analyzing time-series data that 
represent nonstationary and nonlinear 
physical phenomena. The method is 
based principally on the concept of em- 
pirical mode decomposition (EMD), ac- 
cording to which any complicated signal 
(as represented by digital samples) can 
be decomposed into a finite number of 


functions, called “intrinsic mode func- 
tions” (IMFs), that admit well-behaved 
Hilbert transforms. The local energies 
and the instantaneous frequencies de- 
rived from the IMFs through Hilbert 
transforms can be used to construct an 
energy-frequency-time distribution, de- 
noted a Hilbert spectrum. The patent 
expands on the description in the cited 
prior article by explaining underlying 
mathematical principles and describing 
details of implementation. The patent 
also describes, as major elements of the 
method, the options of (1) filtering the 
original signal by combining a subset of 


IMFs and (2) fitting a curve to the fil- 
tered signal — something that it may 
not be possible to do with the original 
signal. 

This work was done by Norden E. Huang 

of Goddard Space Flight Center. Fur- 
ther information is contained in a TSP (see 
page 1). 

This invention is owned by NASA, and a 
patent application has been filed. Inquiries con- 
cerning nonexclusive or exclusive license for its 
commercial development should be addressed to 
the Patent Counsel, Goddard Space Flight Cen- 
ter; (301) 286-7351. Refer to GSC-13817-3. 


Quad-Tree Visual-Calculus Analysis of Satellite Coverage 

The computational burden is less than in a pixel representation. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


An improved method of analysis of 
coverage of areas of the Earth by a con- 
stellation of radio-communication or 
scientific-observation satellites has been 
developed. This method is intended to 
supplant an older method in which the 
global-coverage-analysis problem is 
solved from a ground-to-satellite per- 
spective. The older method is suitable 
for coarse-grained analysis of coverage 
of a constellation of a few satellites, but 
the algorithms of the older method are 
too slow and cumbersome for the large 
scope of the problem of analysis of cov- 
erage of a modern constellation of 
many satellites intended to provide 
global coverage all the time. In contrast, 
the present method provides for rapid 
and efficient analysis. This method is 
derived from a satellite-to-ground per- 
spective and involves a unique combina- 
tion of two techniques for multiresolu- 
tion representation of map features on 
the surface of a sphere. 

The first of the two techniques, called 
“visual calculus,” is one that embodies 
the satellite-to-ground perspective and 
assists in the visualization of global cov- 
erage. In visual calculus, a satellite can 
be regarded as “painting” its field of view 
or other coverage field onto the ground 



Using a Multiresolution Map to analyze a circular region, it is not necessary to resort to full resolution 
everywhere. Large areas with the same value can be represented by a few large squares; smaller 
squares are needed only to resolve details of the boundary. 
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to generate simple maps. These maps 
are two-dimensional projections of the 
globe, with ground tracks and station 
masks for each satellite in a constella- 
tion. These maps are processed further 
to generate composite maps that answer 
a variety of simple and complex ques- 
tions regarding coverage, both for spe- 
cific ground points and for regions. 

For the purpose of computational im- 
plementation, a map is represented by 
an integer function of latitude and lon- 
gitude, analogous to the more familiar 
map color as a function of latitude and 
longitude. More specifically, a map ob- 
ject M is defined by M(l,X):(l,X)^i, 
where l is latitude, X is longitude, and i is 
the integer value at the point in ques- 
tion. The processing of maps to gener- 
ate composite maps and answer ques- 
tions about coverage involves, among 
other things, the use of an algebra in 
which map objects are manipulated by 
unary, binary, Boolean, and other opera- 
tors. Examples of operators include one 
that finds the distance from a given 
point to a map boundary and one that 


finds what percentage of a map has a 
given integer value. 

The second technique essential to the 
present method is the use of a quad-tree 
data structure in implementing the visual 
calculus. For this purpose, the map func- 
tion is redefined to provide not only the 
integer value for a given location but also 
information on how far one must travel 
from that location to encounter a change 
in the map. A point quad tree is well suited 
to a multiresolution representation of 
map features on the surface of a sphere. 

tn essence, a quad tree is a data-com- 
pression construct that takes advantage 
of the facts that (1) a typical map con- 
tains only a few regions with a few values 
and (2) at each point, most values in the 
neighborhood are the same value. In a 
quad tree, squares of pixels that have the 
same value are represented as one large 
square. It does not become necessary to 
do calculations on a subdivision of a 
square until a map boundary subdivides 
the square. Thus, it is not necessary to 
fix resolution at the beginning of an 
analysis. If higher resolution is needed 


for a particular map, it can be achieved 
through finer subdivisions. 

The main reason for choosing a quad- 
tree data structure instead of a pixel- 
based data structure is that the point- 
quad-tree structure requires less 
data-storage space. This computational 
advantage arises because the number of 
subdivisions needed in a given case is ap- 
proximately proportional to the length of 
the boundary between regions, instead of 
proportional to the area of the regions as 
in a pixel representation. It is necessary 
to subdivide only those squares that strad- 
dle a boundary (see figure), and because 
a boundary is one-dimensional, the num- 
ber of subdivisions for the required reso- 
lution scales linearly, rather than qua- 
dratically as when subdividing an area 
with pixels. 

This work was done by Martin W. Lo, 
George Hockney, and Bruce Kwan of Caltech 

for NASA’s Jet Propulsion Laboratory. 

Further information is contained in a TSP 
(see page 1 ). 
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